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INTRODUCTION 


NSWC  has  developed  two  Navstar  Geodetic  Receiver  Systems  (NGRSs)  that  utilize 
Stanford  Telecommunication  Inc.  (STI)  GPS  receivers.  The  first  NGRS  was  operational  in 
February  1979,  and  data  was  gathered  at  several  sites  in  the  United  States  during  that  year. 
Throughout  1980,  the  first  system  was  progressively  modified  so  that  it  and  the  second  system 
(developed  in  1980)  had  as  much  common  hardware  and  software  as  practical.  The  primary 
differences  between  the  two  systems  are  related  to  the  differences  in  the  two  STI  receivers 
employed:  the  first  NGRS  utilizes  a  STI  model  5007  receiver,  while  NGRS-2  utilizes  a  model 
5010.  The  later  model  5010  has  the  capability  of  Lj  pseudorange  observations  and  also 
provides  dual-frequency  6-s  digital  Doppler  count  output,  while  the  5007  does  not. 

Each  system,  which  is  composed  of  the  components  listed  in  Table  1,  is  mounted  in 
three  5-ft  racks.  The  microprocessor  software  and  the  processing  software  that  reads  the 
recorded  data  tapes  were  developed  at  NSWC.  A  more  detailed  system  description  is  given 
in  Reference  1. 


Table  1.  NGRS 

Microprocessor  Controller 

Frequency  Standard 

GPS  Receiver 

GPS  Test  Transmitter 

Time-Interval  Counters 

Full-Cycle  Doppler  Counter 

Digital  Clock 

Antenna 

Tape  Recorder 

Keyboard/Display 

Power  Buffer 


Components 

INTEL  80/10 
HP  5061A-004 
STI  5007  or  STI  5010 
STI  5001A 

HP  5370  and  HP  5328 

NSWC 

HP  5 9309 A 

CKU  CA3207 

KENNEDY  9800 

ADM-3 

ELGAR  UPS  252-1 
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NGRS  DATA  DESCRIPTION 


The  original  intent  of  the  NGRS  was  to  utilize  the  pseudorange  solely  as  a  means  to  obtain 
the  local  clock  correction.  Since  microsecond  accuracy  is  adequate  for  that  purpose,  only  single¬ 
frequency  L!  pseudorange  capability  with  limited  resolution  was  provided  by  NGRS-1 . 

The  single-frequency  pseudorange  measurement  was  upgraded  in  NGRS-1  (1980)  by  replacing 
the  10-ns  time-interval  counter  with  a  20-ps  counter.  This  improved  the  RMS  noise  level  from 
1.5  m  (truncation  limited)  to  70  cm  (receiver  limited).  NGRS-2,  which  utilizes  the  model  5010 
receiver,  has  two  200-ps  counters.  In  all  cases,  the  observation  is  the  time  interval  between  the 
local  6-s  epoch  and  the  received  epoch. 

In  order  to  obtain  ionospheric  corrected  range  difference  observations,  dual-frequency  analog 
Doppler  signals  are  used.  Since  these  signals  are  analog,  a  counter  external  to  the  receiver  is  re¬ 
quired.  A  60-s  nominal  count  interval  was  chosen.  This  method  of  providing  continuous  Doppler 
observations  is  used  in  both  systems.  At  each  local  minute  epoch,  the  current  reading  of  the  L, 
or  L2  Doppler  counter  is  strobed  and  the  time  interval  to  the  next  positive  zero  crossing  is 
measured.  These  two  quantities  are  used  to  construct  the  cycle  count  during  the  previous  1-min 
interval.  The  range  difference  can  then  be  computed  using  the  two  frequency  observations*. 

The  time-of-arrival  difference  for  signals  emanating  from  a  satellite  and  received  at  two 
stationary  (in  the  earth  fixed  frame)  sites  can  be  used  to  determine  the  baseline  connecting  the 
two  sites.  The  data  class  described  is  a  range  difference  in  space  at  a  common  time.  The  NGRS 
are  primarily  designed  to  obtain  range  differences  in  time  at  a  fixed  site;  consequently,  the 
Doppler  count  data  cannot  be  used  directly  in  the  time-of-arrival  difference  problem.  If  the  two 
sites  are  considered  to  be  independent  and  do  not  communicate  in  real  time,  the  most  straight¬ 
forward  arrangement  with  the  NGRS  is  to  plan  the  observations  so  that  the  two  receivers  track  the 
same  satellites  simultaneously  and  record  pseudoranges  and  Doppler  counts  in  the  usual  manner. 
When  the  data  are  later  processed,  the  pseudoranges  from  each  satellite-to-site  combination  can  be 
computed,  evaluated  at  a  common  time,  and  differenced  to  obtain  the  desired  range  difference  be¬ 
tween  the  sites.  Both  the  pseudorange  and  Doppler  data  are  processed  together.  However,  since 
the  Doppler  carries  no  absolute  range  information,  the  evaluation  of  the  smoothed  range  is  de¬ 
termined  primarily  by  the  quality  of  the  pseudorange  obsevations.  The  single-frequency  measure¬ 
ments  are  adequate  for  relative  positioning  if  the  sites  are  not  far  apart,  since  then  both  receivers 
will  see  substantially  the  same  ionospheric  and  tropospheric  refraction  effects. 


FORMULATION 


The  first  adjustment  that  must  be  made  to  the  two  data  sets  is  to  convert  to  a  common  time 
system.  In  this  application,  the  broadcast  ephemeris,  the  broadcast  satellite  clock  corrections,  and 
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the  observed  pseudoranges  can  be  used  to  determine  a  set  of  polynomial  coefficients  that  will 
convert  local  time  to  GPS  time* . 

The  6-s  pseudoranges  are  too  frequent  and  too  noisy  to  be  used  individually,  so  a  smoothing 
polynomial  is  employed  that  also  incorporates  the  range  difference  data.  Only  L,  data  are  used, 
because  NGRS-1  lacks  L3  pseudorange  and  the  baseline  involved  is  very  short.  The  ionospheric 
correction  to  the  L,  range  difference  (derived  from  the  Doppler)  is  computed  and  applied  twice 
to  account  for  the  difference  in  sign  between  phase  (Doppler)  and  group  (range)  delays  introduced 
by  the  ionosphere3.  The  time  tags  attached  to  the  pseudoranges  are  the  local  times  of  reception  of 
the  satellite  epoch ,  while  the  range  differences  are  tagged  at  the  local  epoch.  These  differ  by  the 
amount  of  the  pseudorange  measurement. 

The  smoothing  range  polynomial  is  defined  as  follows: 

r(tk)  =  ao+E  (1) 

i=  1 

In  the  case  of  range  difference,  it  is  r  (t)  -  r(t_  t  ). 

r(tj)~  r(tj_j)=  22  -  toy]  (2) 

i=  l 

These  two  observation  equations  are  used  together  in  one  least  squares  solution  to  solve  for  the 
aj’s.  In  a  typical  10-min  fit,  there  are  100  pseudorange  observations  and  10  range  difference 
observations.  However,  the  range  differences  are  more  accurate  and  are  weighted  by  a  factor  of 
20  to  1  over  the  pseudoranges. 

Similar  10-min  data  spans  from  each  site  are  smoothed  and  the  resulting  coefficients  are  used 
in  Equation  1  to  compute  a  range  at  a  time  near  the  middle  of  the  span.  The  ranges  are  then 
differenced  to  obtain  the  equivalent  time-of-arrival  difference  of  the  satellite  signal  at  the  two  sites. 
Instead  of  fitting  first  and  then  differencing,  the  data  might  have  been  differenced  first  and  then 
fitted.  This  has  the  following  advantages:  allows  the  use  of  a  lower-order  polynomial,  introduces 
less  numerical  error,  and  utilizes  a  single  fit  rather  than  two  fits.  However,  because  the  local  clocks 
were  not  necessarily  synchronized,  the  data  were  not  obtained  simultaneously;  consequently, 
this  option  was  abandoned. 


The  geometry  of  the  relative  positioning  problem  is  illustrated  in  Figure  1.  The  pertinent 
vectors  are  defined  below. 

tT  time  of  transmission 

tR  time  of  reception 

r4(tT )  =  xg  x  +  ys  y  +  zs  z  satellite  position 

rg  (tR )  =  xg  x  +  yg  y  +  zg  z  reference  site  position 

fh(tR )  =  xh  x  +  yh  y +  zh  z  secondary  site  position 

5(tR )  =  rh  (tR )  -  rg  (tR  )  baseline  vector 

r  =  |F(tT  )  -  rg(tR  )j  slant  range,  reference  site 

r  +  C  =  |  r s  (tT  )  -  rh  (tR  )|  slant  range ,  secondary  site 


The  x,  y,  z  coordinate  system  in  Figure  1  is  an  earth  fixed  frame;  therefore,  the  vectors 
F ,  r"h,  and  b  are  independent  of  time.  The  satellite  vector  r^(tT )  is  the  only  time-variable  parameter. 
The  differences  in  range  from  the  satellite  to  the  two  sites  is  represented  by  £  and  is  part  of  the 
observation  obtained  from  differencing  the  two  smoothing  polynomials.  The  remaining  portion 
of  the  observation  is  due  to  differences  in  receiver  biases  and  will  be  represented  by 
Ci  +  Cj(tR  -  tQ).  Thus,  the  observation  equation  is 


L("R)*£  +  £,  +C,(tR  -  to) 


(3) 


The  geometric  part,  C,  is  a  function  o i  the  unknown  baseline  vector  b.  In  the  triangle  GHS, 

(r  +  £)2  =  r2  +  b2  -  2  rb  cos  y  (4) 

The  term  rb  cos  y  is  the  dot  product  b  •  (rf  -  rg ). 

rb  cos  y  =  b  •  (r,  -  rg)  =  xb(x>  -  xg)  +  yb(ys  -  yg)  +  zb  (zs  -  zg)  (5) 

Equation  4  can  now  be  solved  for  C  by  taking  the  square  root  and  substituting  Equation  5. 

C  =  [r2  +  b2  -  2"b  •  (js  -  rf )]%  -  r  (6) 


The  parameters  that  are  to  be  included  in  the  solution  are 

xb’yb’zb’®i>and  22> 

which  represent  the  baseline  vector  and  the  receiver  differences. 

For  each  parameter,  there  is  an  initial  estimate  that  is  to  be  improved  by  using  the  obser¬ 
vations.  The  computed  value  Lc  is  obtained  from  Equation  3,  using  the  initial  estimates.  The 
observed  value  is  LQ .  These  will  differ  by  an  amount  AL  such  that 

L  =  L  +  AL 

O  C 

The  problem  is  linearized  if  Lc  and  AL  are  taken  to  be  the  first  and  second  terms  of  a  Taylor’s 
series  expanded  about  the  initial  estimate.  The  third  and  higher  terms  are  omitted.  The  expanded 
second  term  is 


3LC  3L  b  Lc  9L  3LC 

AL  =  - —  Ax.  +- —  Ayh  +  —  Az.  +  —  A£t  +  —  A£2 

3xb  ”  3yb  3zb  3e>  3C* 


(8) 
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The  required  partial  derivatives  are  obtained  from  Equations  3  and  6: 


3L, 

9xb 


[r2  +  b2  -  2b  •  (r“  -  Fg)]“»  (xb  -  x,  +  xg) 


3L 

—  =  [r2+b2-2b-(Fs-Tg)]-*  (yb-y,  +  yg) 


3L 

ZT  =  [r2+b2-2b-(r|-Fg)]-^(zb-zt  +  zg) 


(9) 


38 1 
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acj 


The  system  of  equations  represented  by  Equation  7  can  be  placed  into  matrix  form  and 
solved  using  weighted  batch  least  squares.  Each  row  represents  one  smoothed  observation. 


<Lo  -  Lc>, 

/3L  \  /3L  \  / 3L  \ 

i  Ix~  /  V  a /  \  IzT /  1  i  ~  4o ) 

Axb 

\0Xb/l  \0yb/l  \azbfl 

<Lo  -  Lc>2 

/ 3L  \  /3L  \  / 3L  \ 

VaT"/  Vav~J  Vaz~/  1  ^2'^ 

Ayb 

\dXb/2  \°Vb/2  \  °Zb/2 

9 

• 

= 

9  9  9  9 

•  99  9 

• 

Azb 

• 

• 

•  99  9 

9  9  9  9 

(Lo  -  LA 

/3L\  /  3L  \  /3L  \ 

(ar)  l  jt)  vst)  1  (t*»'t«>) 

A«s 

'0Xb/n  '  “yb '  n  \0Zb /n 

The  matrix  of  partials  is  denoted  by  A,  and  the  parameters  to  be  improved  are  denoted  by  X. 
The  above  then  reduces  to  the  form 


L  - L  - AX 

o  c 
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Then, 


AtWAX  =  AtW(Lo  -  Lc) 


and 

X  =  (AtWA)_1  AtW(Ld  -  Lc) 


The  matrix  (ATWA)-1  is  the  covariance  matrix. 

The  diagonal  weight  matrix  W  is  introduced  to  give  smoothed  observations  with  low  noise 
more  weight  in  the  solution  than  those  with  high  noise.  The  relative  weights  are  obtained  from 
the  RSS  of  the  RMS  computed  from  the  two  smoothing  polynomials  that  comprise  the  particular 
observation. 

The  solution  vector  X  is  then  added  to  the  initial  estimate  of 

xb ,  yb ,  zfe ,  ii, ,  and  i2 

to  obtain  the  improved  values.  The  standard  deviation  in  this  result  is  obtained  from  the  square 
root  of  the  diagonals  of  the  covariance  matrix.  The  confidence  that  one  places  in  these  standard 
deviations  depends  upon  how  well  the  fluctuations  in  the  observations  agree  with  the  statistical 
assumptions  inherent  in  the  least-squares  process.  Unmodeled  effects  (e.g.,  systematic  errors  in 
Fj  or  variable  instrumental  biases)  can  lead  to  solutions  whose  true  error  is  much  larger  than  the 
computed  standard  deviation.  When  unmodeled  effects  are  known  to  be  present  in  the  data,  the 
standard  deviations  can  only  be  viewed  as  a  guide  indicating  the  geometric  strength  of  the  obser¬ 
vations  and  subsequent  solution. 


OBSERVATIONS 


January  1981,  both  NGRS-1  and  NGRS-2  were  operated  side  by  side  at  NSWC.  Pseudorange 
and  Doppler  data  were  recorded  in  the  usual  manner  while  the  receivers  simulaneously  observed 
the  same  GPS  satellites.  Two  antenna  sites  were  in  use:  a  surveyed  site  and  a  site  approximately 
2  m  south.  The  receivers  spent  time  connected  to  a  common  antenna  and  to  their  individual 
antenna.  Also,  a  common  frequency  standard  was  used,  as  well  as  individual  frequency  standards. 
Table  2  lists  the  number  of  20-min  segments  obtained  from  each  satellite  on  each  day.  For  each 
segment,  a  residual  difference  in  range  is  computed..  The  mean  and  standard  deviations  of  these 
residuals  for  all  segments  on  a  given  day  are  presented.  The  significance  of  these  residuals  is 
discussed  in  the  next  section. 
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Table  2.  Segment  Statistics  For  Each  Days  Observations 


Day  (January  1981 ) 

10 

12 

15 

16 

17 

18 

19 

Antenna 

Common 

Separate 

Common 

Frequency  Standard 

j  Common 

Common 

Separate 

SV-5 

Number  of  Segments 

23 

32 

25 

22 

27 

28 

14 

Mean  (m) 

0.129 

0.138 

0.207 

0.208 

0.112 

0.116 

0.107 

Standard  Deviation  (m) 

0.019 

0.026 

0.038 

0.032 

0.021 

0.024 

0.014 

SV-6 

Number  of  Segments 

14 

8 

11 

14 

14 

14 

13 

Mean  (m) 

0.110 

0.114 

0.218 

0.221 

0.086 

0.091 

0.095 

Standard  Deviation  (m) 

0.009 

0.018 

0.025 

0.042 

0.018 

0.011 

0.015 

SV-8 

Number  of  Segments 

10 

0 

8 

12 

15 

15 

16 

Mean  (m) 

0.109 

0.188 

0.220 

0.090 

0.090 

0.094 

Standard  Deviation  (m) 

0.013 

0.007 

0.025 

0.024 

0.019 

0.012 

SV-9 

Number  of  Segments 

18 

18 

17 

7 

18 

16 

5 

Mean  (m) 

0.158 

0.154 

0.212 

0.265 

0.133 

0.144 

0.149 

Standard  Deviation  (m) 

0.038 

0.032 

0.025 

0.062 

0.031 

0.027 

0.031 

NOISE  PROBLEM 


Unfortunately,  a  problem  with  the  5010  receiver  in  NGRS-2  caused  the  pseudorange 
observations  to  have  a  higher  noise  level  than  comparable  observations  obtained  from  NGRS-1. 
This  fault,  which  was  later  identified  and  corrected,  degraded  the  NGRS-2  pseudorange  data 
somewhat.  The  nature  of  the  degradation  is  illustrated  by  the  Figures  2  through  5. 

Twenty-min  segments  of  data  were  fit  by  a  fourth-degree  polynomial,  as  previously  described. 
Both  L,  6-s  pseudoranges  and  L,  60-s  integrated  Doppler  range  differences  contributed  to  the  fit. 
Figure  2  shows  plots  of  the  fit  residuals  of  two  consecutive  20-min  segments  of  data  from  NGRS-1 . 
Figure  3  shows  similar  residuals  over  the  same  span  from  NGRS-2.  The  diamond  symbols  represent 
range  difference  residuals  from  Doppler,  added  consecutively  to  produce  a  biased  range.  For 
plotting  purposes,  the  bias  was  supplied  by  evaluating  the  polynomial  at  the  time  of  the  beginning 
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Figure  3.  NGRS-2  Residuals  From  20-Min  Smoothing  Polynomial, 
NSWC/Dahlgren,  Day  10, 1981 
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Figure  3.  NGRS-2  Residuals  From  20-Min  Smoothing  Polynomial, 
NSWC/Dahlgren,  Day  tO,  1981  (Continued) 
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Figure  4.  Difference*  in  Pseudorange  (20-Min),  NGRS-2-NGRS-1, 
NSWC/Dahlgren,  Day  10, 1981 
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Figure  5.  Differences  in  Range  Differences,  NGRS-2-NGRS-1, 
NSWC/Dahlgren,  Day  10, 1981  (Continued) 
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of  the  first  Doppler  observation.  The  plus  symbols  are  the  pseudorange  residuals.  In  the  case  of 
NGRS-1,  the  biased  ranges  and  pseudorange  are  consistent,  since  the  residuals  are  centered  about 
zero.  In  the  case  of  NGRS-2,  the  consistency  is  not  as  good. 

In  Figure  3(a)  the  biased  range  residuals  deviate  below  the  zero  axis;  and  in  Figure  3(b), 
they  deviate  above  the  axis.  The  problem  appears  to  be  a  random  short-term  drift  in  the  NGRS-2 
pseudoranges.  Evidence  for  this  has  been  gathered  by  plotting  the  differences  of  the  pseudoranges 
obtained  from  the  two  receivers.  The  differences  should  be  constant  over  the  interval,  because 
both  receivers  were  connected  to  the  same  antenna  and  both  were  being  driven  by  the  same 
frequency  standard.  The  difference  plots  (Figures  4)  show  several  meters  of  scatter  and  a  tilt 
with  positive  slope  (Figure  4a)  and  negative  slope  (Figure  4b).  The  differences  in  range  differences 
between  the  two  receivers  (Figure  5)  show  no  similar  trends;  thus,  the  fault  appears  to  lie  with  the 
pseudorange  bias  from  NGRS-2.  TtK  computed  RMS  for  all  these  residual  plots  is  listed  in  Table  3. 


Tatafr  3. 

RMS  of  Residual  Plots 

Figure 

NGRS 

SPAN 

RMS  (m) 

Smoothing  Polynomii' 

2a 

1 

1 

0.59 

2b 

1 

2 

0.57 

3a 

■> 

1 

0.85 

3b 

-> 

2 

0.88 

Pseudorange  Differences 

4a 

2-1 

1 

0.95 

4b 

2-1 

2 

1.09 

Differences  of 

5a 

2-1 

1 

0.68 

Range  Differences 

5b 

2-1 

2 

0.68 

SOLUTIONS 


DAYS  10  AND  12 

Data  from  the  two  days  in  which  the  receivers  used  common  antennas  and  common  fre¬ 
quency  standards  have  been  processed  in  the  manner  described  in  the  previous  sections.  Figure  6 
shows  the  segments  of  the  satellite  ground  tracks  observed  on  day  10.  The  tracks  for  day  12 
were  similar,  except  that  the  observations  of  SV-8  were  omitted. 
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Figure  6.  Azimuth  and  Elevation  Tracks  of  Satellites 
From  Site  on  Day  10 


All  observations  were  smoothed  over  a  period  of  10  min;  therefore,  they  consist  of  100 
pseudoranges  and  10  range  differences.  The  number  of  10-min  segments  available  on  each  day  is 
listed  in  Table  2.  The  baseline  vector  b  was  purposely  offset  by  4  m  in  each-earth  fixed  compon¬ 
ent  in  these  solutions  in  order  to  demonstrate  that  the  data  produces  corrections  to  b  in  the 
correct  direction.  These  corrections  should  be  -4  m  in  each  component.  The  results  obtained  are 
summarized  in  Table  4.  A  second  solution  was  performed  in  which  the  parameter  fi2  was  con¬ 
strained,  by  a  priori  sigmas,  to  zero.  This  procedure  improved  the  result  on  day  10,  but  worsened 
the  result  on  day  12.  Table  5  presents  this  solution. 

An  additional  experiment  using  just  six  10-min  segments  was  performed.  These  segments 
(Table  6)  were  selected  so  that  they  were  spread  in  azimuth  around  the  site  and  were  nearly 
the  same  on  day  12  as  on  day  10.  Observations  from  SV-8  were  omitted  from  both.  The  solutions 
obtained  are  quite  reasonable  (Table  7)  and  illustrate  the  importance  of  choosing  observations 
that  cover  the  site  from  all  sides.  Similar  solutions  using  data  from  just  one  satellite  (SV-6,  SV-8, 
or  SV-9)  were  poor  due  to  the  absence  of  data  from  two  or  more  quadrants. 
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DAYS  IS  AND  16 


Data  from  the  two  days  in  which  the  receivers  used  common  frequency  standards  but 
separate  antennas  were  processed;  the  baseline  was  to  be  assumed  zero.  In  fact,  the  antenna  for 
NGRS-2  was  placed  2  m  south  and  0.2  m  east  of  the  NGRS-1  antenna.  On  day  15,  the  satellites 
observed  were  SV-5,  6,  and  9;  on  day  16,  SV-5,  6,  and  8  were  observed.  A  solution  using  these 
data  is  presented  in  Table  8.  Separate  solutions  using  only  the  data  from  each  day  are  given  in 
Table  9.  All  of  these  solutions  are  rotated  from  earth-fixed  components  to  a  local  vertical  frame 
with  components  that  are  labeled  North,  East,  and  Vertical,  so  that  comparison  with  the  known 
baseline  is  possible. 

Even  though  the  standard  deviations  indicate  that  the  vertical  component  is  the  best-de¬ 
termined  of  the  three  components,  it  contributes  the  largest  error.  The  excessive  receiver  noise 
during  these  two  days  (as  noted  in  Table  2)  may  have  been  a  contributing  factor  to  the  poor 
results  obtained  in  this  component. 


DAYS  17,  18,  AND  19 

During  the  last  three  days  of  this  experiment,  the  two  systems  operated  from  independent 
frequency  standards,  but  utilized  a  common  antenna.  The  solutions  obtained  were  comparable 
with  those  on  days  10  and  12,  when  a  common  frequency  standard  was  used.  The  values  for  the 
bias  C,  and  drift  fi2  are  larger  than  they  had  previously  been  because  of  the  addition  of  the  relative 
epoch  offset  that  exists  between  the  two  frequency  standards. 

A  combined  solution  using  all  the  data  collected  during  these  three  days  (a  total  of  195 
10-min  segments)  is  presented  in  Table  10.  Again,  a  4-m  offset  in  each  component  was  introduced. 
Separate  solutions  for  days  18  and  19  were  also  made  (Table  11).  Though  the  standard  deviations 
are  smaller  in  the  three-day  solution,  the  magnitude  of  the  error  is  greater  than  that  found  in  the 
single-day  solutions.  This  result  tends  to  indicate  that  the  bias  model  fits  the  data  better  over  a 
short  time  span.  A  result  similar  to  this  has  been  predicted  by  Fell4 . 


* 


f 
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Table  4.  Solutions  Using  All  Available  Data,  Days  10  and  12 


Day 

Component 

Solution 

Sigma 

Solution  -  True 

Units 

10 

X 

-3.570 

0.052 

0.430 

m 

y 

-5.171 

0.109 

-1.171 

m 

z 

-3.448 

0.090 

0.552 

m 

fi, 

38.580 

0.093 

— 

m 

fi2 

17.658  X  10'6 

0.996  X  10"6 

- 

m/s 

12 

X 

-3.955 

0.060 

0.045 

m 

y 

-4.432 

0.130 

-0.432 

m 

z 

-3.830 

0.101 

0.170 

m 

fil 

38.499 

0.105 

— 

m 

fi2 

-1.714  X  10~6 

0.935  X  10" 6 

- 

m/s 

Day  10 

Day  12 

Magnitude  (Solution  -  True) 

1.364 

0.466 

Table  5.  Solutions  Using  All  Available  Data.  Parameter  fi2  Constrained,  Days  10  and  12 

Day 

Component 

Solution 

Sigma 

Solution  -  True 

Units 

10 

X 

-3.543 

0.052 

0.457 

m 

y 

-4.181 

0.094 

-0.181 

m 

z 

-4.279 

0.076 

-0.279 

m 

fil 

38.419 

0.093 

— 

m 

e2 

0.0 

0.0 

- 

m 

12 

X 

-3.978 

0.059 

0.022 

m 

y 

-4.555 

0.112 

-0.555 

m 

z 

-3.719 

0.080 

0.281 

m 

fil 

38.540 

0.103 

- 

m 

fij 

0.0 

0.0 

- 

Magnitude  (Solution  -  True) 


Day  10  Day  12 

0.565  0.622 
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Table  6.  Azimuth  And  Elevation  of  Six  Segments  Selected  From 
Days  10  and  12 


SV-5 


SV-6 


SV-9 


Azimuth  (Degrees) 


Day  10 

Day  12 

310.8 

313.5 

116.9 

115.9 

230.4 

226.5 

211.5 

210.6 

85.4 

85.4 

33.9 

33.9 

Elevation  (Degrees) 


Day  10 

Day  12 

46.2 

47.0 

40.1 

40.8 

61.6 

59.4 

42.6 

41.0 

39.0 

39.1 

21.9 

21.9 

Table  7.  Solutions  Using  Six  Selected  Segments,  Days  10  and  12 


Day 

Component 

Solution 

Sigma 

Solution  -  True 

Units 

10 

X 

-3.805 

0.197 

-0.195 

m 

y 

-3.768 

0.380 

-0.232 

m 

z 

-4.552 

0.376 

0.552 

m 

fit 

37.356 

0.361 

— 

m 

18.736  X  10-6 

3.426  X10~6 

- 

m/s 

12 

X 

-3.557 

0.622 

-0.443 

m 

y 

-5.052 

0.373 

1.052 

m 

z 

-3.192 

0.414 

-0.808 

m 

fit 

39.264 

0.557 

— 

m 

-5.178  X  10“6 

10.559  X10~6 

- 

m/s 

Day  10 

Day  12 

Magnitude  (Solution  -  True) 

0.630 

1.399 

Table  8.  Combined  Solution  Using  All  Available  Data,  Days  15  and  16 


Magnitude  (Solution  -  True) 


Days  IS  and  16 


1.695 


Day 

Component 

Solution 

Sigma 

Solution  -  True 

Units 

15 

North 

-2.122 

0.157 

-.122 

m 

and  16 

East 

0.553 

0.100 

.353 

m 

Vertical 

-1.353 

0.035 

-1.653 

m 

*i 

49.276 

0.133 

— 

m 

C2 

-1.104  X  10"6 

0.557  X  10'6 

— 

m/s 

Table  9.  Separate  Solutions,  Days  15  and  16 


Day 

Component 

Solution* 

Sigma 

Solution  -  True 

Units 

15 

North 

-2.672 

0.507 

-.672 

m 

East 

-0.018 

0.233 

-.220 

m 

Vertical 

1.200 

0.232 

1.200 

m 

50.228 

0.185 

— 

m 

Ci 

-38.980  X  10"6 

12.878  X  10"' 

5  _ 

m/s 

16 

North 

- 1 .844 

0.307 

0.156 

m 

East 

0.869 

0.178 

0.669 

m 

Vertical 

-3.727 

0.104 

-3.727 

m 

c, 

47.254 

0.538 

— 

m 

Ci 

0.366  X  10-6 

6.396  X  10-' 

5  _ 

m/s 

Day  15 

Day  16 

Magnitude  (Solution  -  True) 

1.393  m 

3.790  m 

•Day  15:  SV-5,6,9 
Day  16:  SV-5.  6,  8 
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Table  10.  Combined  Solution  Using  All  Available  Data,  Days  17, 18,  and  19 


Day 

Component 

Solution 

Sigma 

Solution  -  True 

Units 

17,  18, 

X 

-4.512 

0.030 

-0.512 

m 

and  19 

y 

-5.548 

0.047 

-1.548 

m 

z 

-3.327 

0.032 

0.673 

m 

8, 

-3177.808 

0.046 

— 

m 

-136.653  X  10'6 

0.102  X  10- 

6  _ 

m/s 

t 

Days  17,  18  and  19 

Magnitude  (Solution  - 

True) 

1 .764  m 

Table  11.  Separate  Daily  Solutions,  Days  18  and  19 

Day 

Component 

Solution 

Sigma 

Solution  -  True 

Units 

18 

X 

-4.869 

0.052 

-0.869 

m 

y 

-4.432 

0.114 

-0.452 

m 

z 

-4.029 

0.057 

Tl 

-3179.838 

0.202 

m 

«2 

-118.217  X  10-6 

2.206  X  10"6 

m/s 

19 

X 

-3.887 

0.058 

0.133 

m 

y 

-5.402 

0.124 

-1.402 

m 

z 

-3.762 

0.062 

0.238 

m 

-3184.151 

0.408 

— 

m 

«2 

-96.008  X  10"6 

2.334  X  10-6 

- 

m/s 

Dav  18 

Dav  19 

Magnitude  (Solution  -  True) 

1 .303  m 

1 .427  m 
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CONCLUSIONS 


Relative  positioning  on  a  short  baseline  using  single-frequency  pseudoranges  and  range 
differences  from  Doppler  is  permissible  if  the  mission  accuracy  requirements  are  set  at  2-m  or 
greater.  The  time  on  site  need  not  be  great  if  enough  satellites  are  in  view  to  satisfy  the  geometric 
requirements  of  the  desired  solution.  This  is  demonstrated  by  the  data  in  Table  6,  for  which  only  1 
hr  (six  10-min  segments)  of  observations  were  used  to  produce  a  solution  with  errors  comparable 
to  those  obtained  from  10  times  as  much  data. 

The  standard  deviations  in  all  cases  were  highly  optimistic  when  compared  with  the  actual 
errors  obtained.  This  indicates  that  the  RMS  of  the  residuals  to  each  segment  is  not  representative 
of  the  true  errors  involved.  There  appears  to  be  something  akin  to  time-variable  channel  biases 
operating  that  cannot  be  absorbed  by  the  8,  and  fi2  parameters.  Some  method  for  on-line 
calibration,  using  locally  generated  stable  test  signals,  is  needed  in  order  to  derive  corrections  for 
this  interchannel  bias  fluctuation. 

Other  methods  of  treating  these  data  are  capable  of  producing  better  results.  In  particular, 
relative  positioning  using  the  Doppler  data  alone  allows  a  more  precise  solution  to  be  obtained. 
Figures  2  and  3  can  be  used  to  indicate  the  superior  quality  of  the  accumulated  Doppler  vs  the 
pseudoranges  for  this  data  set.  Processing  of  this  same  data  span  using  Doppler  only  has  been 
reported  elsewhere.5 
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